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Abstract 

We present a unified model of compound semiconductor growth based on kinetic Monte Carlo sim- 
ulations in tandem with new experimental results that can describe and predict the mechanisms for the 
formation of various types of nanostructures observed during droplet epitaxy. The crucial features of the 
model include the explicit and independent representation of atoms with different species and the ability 
to treat solid and liquid phases independently Using this model, we examine nanostructural evolution 
in droplet epitaxy The model faithfully captures several of the experimentally observed structures, in- 
cluding compact islands and nanorings. Moreover, simulations show the presence of Ga/GaAs core-shell 
structures that we validate experimentally A fully analytical model of droplet epitaxy that explains the 
relationship between growth conditions and the resulting nanostructures is presented, yielding key insight 
into the mechanisms of droplet epitaxy 

1 Introduction 

Nanotechnology involves the manipulation of matter at the nanoscale in order to take advantage of the 
physical properties of materials which arise by a fine tuning of shape and size. Due to quantum con- 
finement effects, three dimensional semiconductor nanostructures can behave as artificial atoms and, like 
natural atoms, show^ a discrete spectrum of energy levels UJ. In contrast to actual atoms, the electronic 
properties of quantum nanostructures can be finely tuned by adjusting structural parameters such as size, 
composition and morphology. The latter parameter is the most relevant for the control of quantum nanos- 
tructure electronic properties, as tiny variations in morphology or composition can cause dramatic changes 
on the band structure [2J. 

One of the most common methods for the fabrication of quantum nanostructures is the growth of lattice- 
mismatched III-V semiconductor materials via the Stranski-Krastanov (SK) mode. SK growth exploits the 
self-assembly of pyramidal-like quantum dots (QDs), driven by the relaxation of lattice mismatched strain 
accumulated in the epilayer. Despite the high success of the technique, which led to fundamental physical 
understanding and a variety of applications, the available design degrees of freedom remain limited. The 



1 



precise engineering of size and shape of QDs via SK self-assembly remains problematic, thus limiting the 
possibilities for the design of arbitrary electronic and optical properties. 

To overcome the SK growth limitations, a fully kinetic limited growth procedure called droplet epitaxy 
(DE) was introduced IHlllI. Unlike SK self-assembly, DE does not rely on strain for the formation of three 
dimensional nanostructures. Instead, DE is based on the sequential deposition of III and V column ele- 
ments at controlled temperatures and fluxes. An initial deposition of metallic group III such as Ga in a 
group V free environment leads to the formation of nanometer size liquid droplets. Subsequent exposure 
of the liquid droplets to the group V vapor causes them to crystallize into quantum dots ||5l[6|. A variety 
of nanostructures may be obtained using this technique, ranging from compact quantum dots, quantum 
rings, and core shell structures ISVfTTl. While experimental observations evoke likely mechanisms, a unified 
model for the formation of various nanostructural morphologies via DE is still needed. Understanding the 
precise mechanisms behind their formation is critical to controlling the various properties such as shape, 
size, and composition, which in turn affect the macroscopic behavior of devices based on such quantum 
nanostructures. 

Kinetic Monte Carlo (KMC) simulation is a useful tool for developing and testing models of epitaxial 
growth and have shed light on the important processes for a wide range of scenarios. Previous KMC simu- 
lations of GaAs systems modeled homoepitaxial growth on a cubic lattice. Among the earliest simulations, 
a single component was used to model epitaxial growth, varying energy parameters to account for the ex- 
pected effect on diffusivity by different Ga/As deposition ratios iTTZll . Two-component simulations have 
also been performed IIT3UT51 to study aspects of film growth such as the transition from two dimensional 
island formation to step-flow growth. Such simulations, however, are necessarily over small areas and are 
valid only for a limited range growth parameters. In most cases these simulations are unable to model 
non-stoichiometric growth; that is, growth when the group V to group III ratio in the growing film is not 
unity, as in the case in droplet epitaxy. 

In this paper, we present experimental and simulation results that shed light on the exact mechanisms 
behind the formation of nanostructures in droplet epitaxy. We introduce a KMC model (Section O that is 
capable of simulating the formation of nanostructures at all stages, from homoepitaxial growth, to liquid 
metal droplet formation and subsequent crystallization. The model represents Ga and As atoms explic- 
itly, and treats liquid and solid phases independently. Although atomistic, the simulations form large-scale 
global structure consistent with experimental results in a reasonable amount of time. As an initial valida- 
tion of the model, GaAs homoepitaxy simulations are compared with analogous experiments, accurately 
reproducing a surface termination phase diagram (Section|3|. We then turn our attention to the main thrust 
of the paper, the nanostructural dependence on growth conditions during droplet epitaxy. We present new 
experimental results exhibiting this dependence, along with corresponding simulations (Section |4]). The 
model can capture the broad range of nanostructures observed in the experiments with the correct qualita- 
tive dependence on growth parameters. Moreover, simulations predict the presence of Ga/ GaAs core-shell 
structures, which are difficult to observe unambiguously in experiments. The existence of both monocrys- 
talline and polycrystalline shells suggest two independent mechanisms of their formation: a morphological 
instability of the crystallization front or nucleation at the vapor-liquid interface. We show by simulation that 
nucleation-induced shell structures may be recrystallized into fully crystallized GaAs islands by annealing 
in high temperature (Section [4.2. In the case of the instability-driven shell formation, we provide sim- 
ulation and experimental evidence to suggest the presence of a Mullins-Sekerka instability (Section [4.2.21) . 
Lastly, we develop a fully analytical model (Section O that describes the existence of the structures ob- 
served and their dependence on growth conditions. The theoretical model agrees well with simulation and 
experimental results. 

2 Kinetic Monte Carlo model 

Prior work in simulating GaAs systems focused on homoepitaxial growth and studied associated phenom- 
ena such as step density [12| and growth modes of GaAs films [13J. In these simulations, surface diffusion 
of adatoms played the central role, and because of the stoichiometric nature of epitaxial growth, a simple 
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cubic lattice and the solid-on-solid constraint sufficed in modeling this process. Droplet epitaxy however, 
poses several issues that cannot be captured by earlier models. During DE, the system is inherently non- 
stoichiometric, so care must be taken to model liquid droplet formation. Processes other than surface diffu- 
sion, such as events within liquid and at the liquid-solid interface, play a key role and cannot be captured 
by a simple solid-on-solid model. Lastly, DE results in large-scale nanostructures on the order of 10s to 100s 
of nm so that the simulations must be performed efficiently within large domains. 

We simulate the homoepitaxial growth and droplet epitaxy experiments with a KMC model that at- 
tempts to address the relevant processes in DE. In the model, atoms occupy positions on a 1+1 dimensional 
analog of the zinc blende lattice illustrated in FigureHl Each position on the lattice is adjacent to four nearest 
neighbors and four next-nearest neighbors. Lattice sites are vacant or occupied by either a Ga or As atom. 
The 1+1 dimensional approach naturally hinders our model to catch some subtlety in quantum nanos- 
tructure shape control that are possible in droplet epitaxy thanks to fine control of surface reconstruction 
dynamics |[T6l . 

In each Monte Carlo step, one of four local events alters at most two positions on the lattice. Atoms may 
desorb from the surface (A in Figure [l]), diffuse on the surface (B), exchange with neighboring atoms (C), 
or adsorb onto a vacant position on the surface (D). The model does not enforce a solid-on-solid constraint. 
Instead, atom configurations are required to be connected, i.e. there must exist a chain of atoms (through 
nearest neighbor bonds) between any two atoms in the configuration. While connectedness is inherently 
a global property, the requirement is approximated by enforcing the property on all local neighborhoods. 
Any event resulting in a disconnected configuration is disallowed. 




Figure 1: KMC Events. Simulations are run on a 1+1 analog of GaAs zinc blende structure. Ga atoms are 
colored red. As atoms are green. The model allows for atom desorption (A), surface diffusion (B), atom- 
atom exchanges (C), and atom deposition (D). The simulations are explicitly multi-species and atomistic, 
and evolve via events on individual Ga or As atoms. 

Event rates are determined by simple nearest and next-nearest neighbor bond-counting. Due to the 
multi-species nature of the model, several bonding energies must be specified. Specifically, nearest and 
next-nearest neighbors are bonded by energies that depend solely on the species of the two connected 
atoms. Nearest neighbor bonding energies are given by three parameters 7 (Ga, Ga) , 7 (Ga, As) and 7 (As, As) . 
Next-nearest neighbor bonds are only assigned between atoms of the same species and are denoted jnn (Ga, Ga) 
and jnn (As, As) . In this way, Ga and As only interact through nearest-neighbor bonds, whereas Ga-Ga and 
As- As interact through nearest and next-nearest neighbor bonds. The parameter values used in the sim- 
ulations are summarized in Table [H These energy values imply that the Ga-As bonds are energetically 
preferred, while weak As-As bonds effectively eliminate excess As by desorption. The bonding energies 
were calibrated from the homoepitaxy (Section O and liquid droplet simulations (Section ID) to match qual- 
itative experimental assumptions and observations. Within the homoepitaxy simulations, the Ga-As bond 
strengths were tuned in order to observe a transition from rough island formation in the low temperature 
regime to a smooth step-flow growth mechanism in high temperatures fV7\. In the case of liquid droplet 
simulations, we tune the Ga-Ga bond strength to match liquid droplet statistics such as droplet width and 
number density. 

Monte Carlo transition rates r(X, Y) from atom configuration X to configuration Y are taken to be of an 
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Table 1: Pairwise nearest and next-nearest neighbor bonding energies used throughout the paper. 
Arrhenius form: 

r(X,y) =]^o^-^«(^'^)/^^^ 

where i^o = 10^^ sec~^ is a constant rate prefactor, T is temperature, and Ea{X,Y) is the activation energy 
for a particular event. 

For surface diffusion, Ea is simply the sum of the bonding energies about the diffusing atom. For the 
desorption of an As atom, the activation energy is the sum of the bonding energies about the desorbing 
atom and an additional desorption barrier = 1.1 eV, its value calibrated by homoepitaxial simulations 
described in Section O The desorption of Ga is disallowed. In the case of atom-atom exchanges, which 
occur within a liquid droplet, we use a base energy barrier for the As diffusion barrier in liquid Ga, 
described in more detail below. To facilitate exchanges occurring at the liquid-solid interface, for which 
we must take solid bonds into account, we define an intermediate state in which the exchanging atoms 
are replaced by an intermediate species H. Bonding energies between Ga or As atoms and H atoms are 
obtained by averaging over H = Ga and H = As, e.g. 

7(Ga,H) = ^(7(Ga,Ga)+7(Ga,As)). 

The activation energy Ea of an atom-atom exchange at the liquid-solid interface is then taken to be the sum 
of the base barrier ejj and the change in bonding energies between the original and the intermediate states. 

For example, we consider the attachment and detachment of an As atom onto the liquid-solid interface, 
as illustrated in Figure |2l We denote the states before and after the detachment as X, Y, respectively, and 
the intermediate state as defined above by J(X, Y). The change in energy £(X) — £( J(X, Y)) between states 
X and J(X, Y), describing the additional barrier for attachment can be computed from a small number of 
relevant bonds, indicated by black lines in the figure. The total activation energy for attachment £^ = 
€d + - £(I(X, Y)) is then given by 

£a = + 47 (Ga, As) + 47 (Ga, Ga) + 2jnn (Ga, Ga) 

-87(Ga,H) -57Hn(Ga,H) -7^^(As,H), 

1 1 
= - ^Tnn (Ga, Ga) - -jnn (As, As), 

= 0.5 eV. 

We may calculate the activation energy for detachment £d = + £(^) ~ £(^(^/ Y)) in a similar manner: 

£0 = ^0+ 47 (Ga, As) + 47 (Ga, Ga) + 3jnn (Ga, Ga) + jnn (As, As) 
-87(Ga,H) -57„n(Ga,H) -7„h(As,H), 

1 1 

= + -7nn(Ga,Ga) + -7nn(As,As), 

= 0.9 eV. 

The forms of the attachment and detachment barriers and E^ indicate what physically occurs during the 
transitions, namely the formation/ removal of one next-nearest neighbor Ga-Ga bond and one next-nearest 
neighbor As- As bond. 
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Figure 2: Attachment and detachment events at the liquid-solid interface. This figure illustrates the attach- 
ment and detachment of an As atom in liquid Ga onto and from a perfectly flat liquid-solid interface, along 
with the intermediate state for the transitions. The intermediate species H is colored blue. The black lines 
indicate the relevant bonds that contribute to the change in energy between the initial and intermediate 
state. 



2.1 Liquid Barriers 

Because we wish to simulate liquid droplet epitaxy, our model identifies liquid neighborhoods in which 
events are assigned different energy barriers to account for the salient physics within a liquid droplet. The 
diffusion of As through a liquid droplet occurs atomistically as an exchange event between the diffusing 
As atom and a Ga atom neighbor. Such an atom-atom exchange event is deemed a liquid event if there is 
at most one As atom total within the neighborhood of the exchanging atoms. 

Further, such exchanges in the liquid are categorized as occurring on either the surface or bulk of a 
droplet. A surface liquid exchange event is one in where there is at least one vacancy within the neighbor- 
hoods of the exchanging atoms. These events represent the incorporation of an As atom from the surface 
into the bulk of the droplet or vice versa and are given a constant activation energy barrier ej = 1.05 eV. The 
particular value of ej is slightly larger than the 1 eV barrier for the diffusion of an As adatom on the surface 
of a liquid droplet. A liquid event that does not occur on the surface is categorized as a liquid bulk event 
and is assigned the energy barrier ejj, which represents the diffusion barrier of an As atom within the bulk 
of the liquid. The value = 0.7 eV was calibrated by the droplet epitaxy simulations described in Section 
m Only after this model was completed did we become aware of the experimental work of Gorokhov et 
al. IITSl who reported the energy barrier for diffusion As in Ga of approximately 0.7 eV. 

2.2 Detailed Balance 

In order to discuss detailed balance, we consider the behavior of our model in the absence of deposition 
and desorption |[T9l . Under such conditions, we assume the incorporation of an As atom from the surface 
of a droplet into the liquid bulk to be a rare event. The rates of the remaining events (i.e. surface diffusion 
and atom-atom exchanges both within the liquid and at the liquid-solid interface) satisfy detailed balance 
with respect to the Boltzmann distribution 

where E{X) is the energy of the state obtained by bond counting, using bonding energies in Table [l] For 
these events X ^ Y, the activation energies are of the form 

Ea = £(x) - £(j(x,y)) + e(x,y), (1) 

where J(X, Y) is a transition state between X and Y and e(X, Y) is an additional barrier. Detailed balance 
follows since we have required that I (X, Y) and e(X, Y) be symmetric with respect to their arguments. 
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For surface diffusion, J(X, Y) is the state with the diffusing atom removed and e(X, Y) = 0. Because 
the reverse event Y ^ X is also a surface diffusion event, it follows that J(X, Y) = J(Y, X) and e(X, Y) = 
= e(Y, X). For atom-atom exchanges, J(X, Y) is the state where the exchanging atoms are replaced 
with intermediate species H, and hence J(X, Y) = J(Y, X), and e(X, Y) = €£,. Note that it was claimed 
above that for an As atom diffusing through liquid Ga, an activation barrier ejj was used without any 
further contribution to the barrier from bond counting. This is consistent with the generalized description 
of the activation energy in equation ([ij, because in this case £(X) = £(J(X, Y)) = £(Y), i.e. no bonds are 
formed/ removed during the diffusion of As in liquid Ga. 

2.3 Implementation 

Because of the time and length scales associated to the nanostructures we wish to simulate, an efficient 
implementation of the KMC algorithm is important. A rejection-free sampling of rates is achieved by the 
standard binary-tree data structure containing partial sums of the rates [20], achieving 0(log N) sampling 
by binary search and 0(log N) tree updates, where N is the number of atoms. As the algorithm proceeds, 
event rates for several atoms must be recomputed at every step. A rate memoization scheme is employed 
to remove repeated rate calculations by taking advantage of the recurrent nature of local neighborhoods 
within KMC [|21||. This procedure leads to significant performance gains. 

3 Homoepitaxy 

As an initial application and validation of the simulations, we model GaAs substrate growth at various 
growth conditions. It is known that the surface reconstruction of a GaAs substrate depends on both the 
temperature and the relative deposition rates of Ga and As ||22|. In 1+1 dimensions, simulations cannot 
reproduce surface reconstructions. Instead, we measured surface termination and its dependence on tem- 
perature and incoming deposition rates. The simulation results were then compared to experimental data. 

Experiments on GaAs(OOl) were carried out using molecular beam epitaxy (MBE) equipped with a solid 
source for Ga and a valved source for AS4. Substrates were initially heated under AS4 overpressure to 
desorb the native oxide layer. A 75 nm GaAs buffer layer was then grown at a substrate temperature 
of 590 °C with a Ga deposition rate of 0.4 monolayers per second (ML/ sec) and an AS4 rate of at least 2 
ML/ sec, as measured by reflection high energy electron diffraction (RHEED) oscillations. After growing 
the GaAs buffer, the substrate temperature was fixed between 460 °C and 610 °C. GaAs was then grown 
with a Ga rate of 0.6 ML/ sec while the As rate was slowly reduced from 3.9 ML/ sec to 1.0 ML/ sec. During 
the growth, the RHEED pattern was monitored for the change from an As-terminated reconstruction to a 
Ga-terminated reconstruction. Consistent with the expected result, at high substrate temperatures the As 
rate required to maintain an As-terminated surface increases with temperature. At lower temperatures, the 
necessary As rate becomes invariant with temperature. 

To simulate GaAs substrate growth, both Ga and As were deposited simultaneously on initially flat. 
As-terminated substrate. The Ga deposition rate f ga was fixed at 0.37 ML/ sec, while the As deposition 
rate f^s was varied so that the deposition ratio fAs/^Ga ranged between .5 and 10. The temperature was 
varied between 427 °C and 727 °C. Five monolayers of total material was deposited, and surface Ga con- 
centration was measured at regular intervals during the deposition of the last two monolayers. Figured is a 
surface termination phase diagram as a function of deposition rate and temperature for both experimental 
and simulation results. Red squares indicate the conditions where simulations show a predominantly Ga- 
terminated surface, while green circles are conditions yielding a predominantly As-terminated surface. Ex- 
perimentally determined transition from the As-terminated 2 x 4 to the Ga-terminated 4x2 reconstructions 
as determined by RHEED are shown as blue diamonds and indicate a good agreement between simulation 
and experimental data. 

The relevant parameter that controls surface termination is the As desorption barrier jij^g. Varying this 
parameter effectively shifts the above phase diagram horizontally. The value ji^s = 1-1 eV is fitted to 
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experiments. This value, combined with the energy values in Tabled] yield a total activation energy of 



2j{Ga, As) + jnn {As, As) + jiAs = 2.2 eV, 

for the desorption of an As adatom from a Ga terminated substrate, comparing favorably to experimental 
results II23II24I . Moreover, the specific parameter value does not significantly impact the qualitative shape of 
the phase diagram. That is, independent of jiAsf the simulations capture a constant critical deposition ratio 
in the low temperature regime and its transition to an increasing critical ratio as temperature increases. 

An analytical expression for the boundary between Ga-terminated and As-terminated substrates can 
be determined. During homoepitaxial growth, several processes occur, including surface diffusion and As 
desorption, but to first order we may approximate the system in a quasi-static deposition/ desorption- 
limited regime. In this regime, the transition between the As and Ga-terminated surface occurs when the 
amount of Ga on the impinging upon the surface (given by Fca) is equal to the net rate of As growth (given 
by f As ~ ^desorb)r where Rdesorb is the desorption rate of As, taken as the harmonic average of the desorption 
rates for an As adatom on a Ga-terminated surface (with desorption energy barrier 2.2 eV) and a As atom 
on a Ga-terminated surface with one nearest or next nearest neighbor (with desorption barrier 2.3 eV). This 
boundary is given by the equation 



^As -1 I -^desorb 



= 1 + 



(2) 



and agrees well with both simulation and experiments. 
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Figure 3: Substrate termination phase diagram as a function of deposition ratio and temperature, obtained 
from simulations and experiments. Those growth conditions resulting in a mostly Ga-terminated substrate 
are in indicated by red squares, while green circles label As-terminated ones. The blue points above indicate 
the conditions where the transition from Ga to As termination occurred experimentally. The blue curve 
indicates the boundary between Ga and As-terminated given by equation |2l 



4 Liquid Droplet Epitaxy 

The main thrust of this work was to accurately capture the processes relevant to liquid droplet formation 
and crystallization. The experiment and simulations proceed in two main steps, which are at the basis 
of the droplet fabrication of GaAs/AlGaAs quantum nanostructures \M- First, Ga is deposited onto the 
substrate by MBE, forming liquid Ga droplets. Next the droplets are crystallized by the introduction of an 
As flux. Depending on growth conditions such as temperature and As overpressure, the actual morphology 
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upon crystallization can vary (SI. The observed nanostructures may range from fully crystallized compact 
quantum dots, to nanorings and even etched holes |[7l[8l [T6ll25ll26ll . 

4.1 Quantum Dots and Nanorings 

4.1.1 Experimental methodology 

In order to experimentally assess the exact dependence of droplet epitaxy nanostructures on the growth 
condition, we prepared a matrix of samples, sampling the typical droplet epitaxy ranges for substrate tem- 
peratures and As4 beam equivalent pressure (BEP). The samples were prepared following the standard 
procedure for droplet epitaxy L6,2Z*i28J. After the oxide removal and the GaAs buffer layer growth (1 ]im 
on thick at 580 °C) to ensure the atomic smoothness of the surface, a 200 nm thick Alo.3Gao.7As barrier layer 
is grown at the same temperature of 580 °C. The presence of an AlGaAs barrier permits the use of the 
fabricated structures, once capped with another AlGaAs layer, as morphology controlled quantum nanos- 
tructures. We do not expect the Al content of the barrier to influence the obtained GaAs nanostructure 
morphology as the AlGaAs barrier is buried under a 1.75 ML GaAs top layer which forms on the surface 
during the following Ga deposition step 1 6.29 1 . 

Prior to the deposition of Ga, the As cell is closed and the background pressure reduced below 10~^ Torr. 
Ga droplets are then formed on the substrate surface by supplying 2.5 ML of Ga at 350 °C with a deposition 
rate of 0.08 ML/ s. Many nearly hemispherical Ga droplets form at a density of 8 x 10^ cm~^. The average 
droplet diameter and height are 50 nm and 20 nm, respectively. The same Ga droplet preparation procedure 
is used for each sample in order to assure the same droplet density before the crystallization procedure. We 
then explored the nanostructure fabrication parameter space by varying the substrate temperature and the 
As4 BEP used for the crystallization of Ga droplets. Nine different samples are prepared by systematically 
varying the temperature between 150 and 350 °C and the AS4 BEP between 5 x 10~^ and 5 x 10~^ Torr. 
The crystallization time is kept constant (10 min) to ensure the complete reaction of the Ga contained in 
the droplets with As atoms. The morphological dependence of the fabricated nanostructures, measured by 
Atomic Force Microscopy (AFM) on temperature and As overpressure is summarized in Figure IH From 
this data, it is clear that compact islands form at low temperature, and as the As overpressure is reduced 
or the temperature is increased, the nanostructures become rings. At the highest temperature and lowest 
As overpressure, only holes remain in the place of the droplet. From large area scans (not shown here) the 
number of GaAs nanostructures per unit area are in excellent agreement with the original droplets density. 
Each droplet was thus transformed into a GaAs nanostructure at the end of the crystallization procedure. 

4.1.2 Simulations 

The simulations were performed in a similar manner to the experiments outlined above. However, in order 
to simplify the analysis, the temperature was maintained constant throughout each simulation. This is 
justified by the observation that the Ga deposition temperature affects droplet density but does not play a 
relevant role in determining droplet crystallization dynamics. Simulations proceed as follows: 

1. Ga atoms are deposited on a flat. As-terminated GaAs substrate at a rate of 0.1 ML/ sec at temperature 
T until 4.0 monolayers are deposited; 

2. The system is then annealed for 60 seconds in the absence of deposition; 

3. After annealing, an As flux is introduced by the deposition of As at a rate of f^g ML/sec until the 
system attains equilibrium. 

The growth parameters T and f^g were varied in order to study their effect on the resulting morphology. 
The temperature T ranged between 150 °C and 350 °C, while the As deposition rate ranged between 0.1 
ML/ sec and 4 ML/ sec. Within this range of growth parameters, we are able to simulate the formation of a 
variety of nanostructures similar to those observed experimentally. 
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Figure 4: Droplet epitaxy experimental results showing typical nanostructures observed over a range of 
substrate temperatures and AS4 BEP. 

During the first phase of the simulation, Ga atoms are deposited on a flat, initially As-terminated GaAs 
substrate. The first monolayer of Ga deposited is consumed in creating a layer of Ga, resulting in a Ga 
terminated substrate. Afterwards, the remaining Ga atoms diffuse along the surface and eventually nucle- 
ate hemispherical droplets as the system attempts to minimize the vapor-liquid interface. During droplet 
formation, simulations show that liquid Ga etches into GaAs, and the amount of etching is regulated by 
temperature, as illustrated in Figure |5l Higher temperatures result in more significant etching. This is in 
agreement with experimental observations II251I3QI . As the droplet etches into the substrate liquid Ga atoms 
displace substrate As atoms, which subsequently attach near the triple junction. In addition, some of the 
displaced substrate material is wicked out of the droplet in a step-flow growth mode. The relevant model 
parameter controlling the effect of etching is the additional barrier for atom-atom exchanges, Cd. The value 
€£, = 0.7 eV was selected to fit qualitative experimental observations on the amount of etching occurring at 
various temperatures. 

Once the liquid droplet has formed and come to equilibrium, an As flux is introduced to initiate crys- 
tallization. Arsenic atoms deposited near or on the droplet diffuse through the liquid quickly IITSl and 
attach most typically near the triple-junction. Such crystallization results in a growing GaAs front and the 
droplet is crystallized inward. If no nucleation occurs at the vapor-liquid interface and the GaAs fronts 
coalesce, a fully crystallized quantum dot forms in place of the liquid droplet. Such is the case for moderate 
temperatures and deposition rate. 

Figure[6]is a sequence of simulation snapshots illustrating the crystallization of a liquid droplet resulting 
in a quantum dot (left panel), along with analogous AFM images of the GaAs fronts obtained experimen- 
tally (right panel). The simulation images in the figure illustrate a typical quantum dot grown at T = 275 °C 
and = 0.06 ML/sec. The general trend is that a Ga drop forms once enough Ga has been deposited 
on the surface (Figure [6K)/ followed by crystallization near the vapor /liquid /solid triple junction upon 
exposure to As flux (B,C). As crystallization progresses, the liquid Ga is consumed, resulting in a fully 
crystallized quantum dot (D). Experimental images (Figure [H right panel) were obtained from individual 
samples prepared according to the experimental procedure outlined above, varying As exposure time. The 
unreacted liquid Ga was removed from the samples before imaging by selective wet etching [8J, thus show- 
ing GaAs fronts at various stages during crystallization. The AFM images confirm the growth mechanism 
observed in the simulations during crystallization. 

Besides compact quantum dots, other nanostructures are predicted by the simulation by considering 
a broad range of temperatures and As fluxes. For example, if F^g is sufficiently low or T is sufficiently 
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A T = 200°C 




Figure 5: Liquid droplets grown at T = 200, 300, 350 °C and Fca = 0.1 ML/sec. Here and throughout the 
paper, Ga and As atoms initially belonging to the substrate are colored red and green, respectively. Ga and 
As atoms deposited throughout the simulation are colored purple and blue, respectively. 

large, the simulations show the formation of nanorings upon crystallization. The exact morphology of the 
nanorings is sensitive to the growth conditions. Figure [71 shows the morphological dependence of the rings 
on f^s, fixing T = 375°C. At low As deposition rate (f^g = 0.1 ML/sec), broad and short nanorings form. 
As F^s is increased the nanorings become more compact and taller so that at f^g = 0.4 ML/ sec the resulting 
structure resembles a "pitted'' quantum dot. This compares well to the experimental results in FigurelH e.g. 
the transition in structure between Figured]: and Figure |4]f as As4 BEP is lowered. Figure Ht resembles the 
pitted quantum dot structure of Figure [Zt- If the BEP is lowered, the resulting nanostructure in Figure St is 
a broad and shallow disk surrounding a pit, resembling Figure [Zh- Similar structures and their dependence 
on both As deposition rate and temperature have been reported in the literature 11011311 . 

4.2 Core-Shell Structures 

In addition to quantum dots and nanorings, simulations show the existence of Ga/GaAs core-shell struc- 
tures. These structures consist of liquid Gallium being completely surrounded by GaAs. In the low tem- 
perature/high deposition rate regime, the GaAs shells are polycrystalline. However, in a higher tempera- 
ture/lower flux regime, the simulations show the formation of a shell in registry with substrate. It will be 
shown below that the first case is the result of nucleation of GaAs at the vapor-liquid interface whereas the 
second case results from a Mullins-Sekerka instability of the crystallization growth front. 

4.2.1 Nucleation 

In the low temperature and high As deposition rate regime nucleation of GaAs clusters near the vapor- 
liquid is significant. This results in the formation of a polycrystalline GaAs shell surrounding a liquid Ga 
core, as illustrated in Figure [U In the figure, a liquid droplet grown at T = 150°C is crystallized by an 
As flux, deposited at a rate of 0.8 ML/ sec. Nucleation at the vapor-liquid interface occurs within seconds 
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Figure 6: Left panel: model snapshots of liquid droplet crystallization at various times for T = 275 °C and 
^As = 0-06 ML/sec resulting in a compact quantum dot. Ga and As atoms from the original substrate 
are colored red and green, respectively. Ga atoms deposited to form liquid droplets are colored purple, 
while As atoms deposited during crystallization are blue. Right panel: AFM images of the GaAs growth 
fronts in partially crystallized droplets after 10 seconds (top), 40 seconds (middle) and 90 seconds (bottom). 
Crystallization was obtained at T = 150°C and 5 x 10~^ Torr As pressure. 



A Fa, = 0.1 ML/sec 




Figure 7: Nanorings formed at T = 375°C and Fas = 0.10, 0.20 and 0.40 ML/sec. 
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upon crystallization (Figure |8p). The liquid core in the final configuration (C) is completely surrounded by 
a GaAs shell after 2.4 seconds, preventing any further crystallization of the liquid. 




Figure 8: Simulation snapshots at times f = 0, 1.8,2.4 seconds after crystallization, T = 150°C and f^g = 0.8 
ML/sec. 

The presence of polycrystalline GaAs, with a high number of grain boundaries and stacking faults that 
this implies, prevents the possibility that such islands could act as efficient emitters and hence are unde- 
sirable. Annealing at higher temperature may remove such defects because grain boundaries and stacking 
faults provide fast diffusion paths for the liquid Ga trapped within the shell. Such paths are accessible 
at higher temperatures. Moreover, thermal fluctuations of the nuclei can effectively serve to dissolve the 
shell, which are characteristically thin in this regime. Therefore such configurations may be annealed at 
a high temperature to remove defects. Indeed, the simulations bear this out. Figure [9] shows a quantum 
dot with polycrystalline shell resulting from the crystallization of a liquid droplet at an As deposition rate 
of f^s = 0.80 ML/sec and temperature T = 150°C . The dot is then annealed at a higher temperature 
T = 350 °C. The initial configuration (|9j\) shows the droplet prior to recrystallization. Temperature is then 
increased, maintaining the same As flux. Ga atoms move along grain boundaries toward the surface, re- 
sulting in a broadening of the dot. The thin shell dissolves, resulting in liquid Ga exposed to As (B, C). 
Within two seconds, the droplet becomes fully crystallized (D) into a shallow GaAs island, absent of any 
defects. 

These simulations show that liquid cores arising due to nucleation at the vapor-liquid interface, which 
would be detrimental to optical and electronic properties of quantum dots, can be eliminated by annealing. 
This expands the parameter space available for the formation of quantum dots formed via droplet epitaxy, 
allowing the fabrication of optical quality dots even at low temperatures. Even if not strictly related to 
the presence of a liquid core, but rather to a low quality of the crystalline nature of the fabricated nanos- 
tructures, post-crystallization in-situ annealing temperatures have been experimentally demonstrated to 
strongly increase the optical quality of droplet epitaxy GaAs/ AlGas dots, reducing individual dot emission 
lines |32|, allowing single photon emission up to liquid nitrogen temperatures f33| and increasing the dot 
emission decay time. 
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Figure 9: Simulation snapshots of a quantum dot annealed at high temperature at time t after tempera- 
ture was increased. (A) The dot after exposure to As deposition at = 0.80 ML/ sec and temperature 
T = 150°C. This results in a polycrystalline GaAs shell trapping a liquid Ga core. (B)-(D) Temperature is 
increased to T = 350 °C and the atoms rearrange in order to fully crystallize the liquid core. 



4.2.2 MuUins Sekerka Instability 

As described above, simulations show that crystallization performed at low temperature and high As flux 
result in nucleation of polycrystalline GaAs shells surrounding a liquid Ga core, resulting from nucleation 
of GaAs at the vapor-liquid interface. In the high temperature regime, such surface nucleation does not 
occur. However, for crystallization at sufficiently high fluxes and high temperature, simulations show the 
existence of liquid Ga core structures surrounded by GaAs shells in registry with the substrate. As a conse- 
quence, such shells are monocrystalline and result from a mechanism separate from surface nucleation. 

By examining the formation of such structures in the simulations, we propose that these shells are driven 
by an instability at the liquid-solid growth front. Simulation snapshots in Figure [lO] illustrate the growth 
mechanism behind this. When this phenomenon occurs no surface nucleation is observed; instead, the 
growth of the GaAs front undergoes an instability at the liquid-solid interface with dendrite-like structures 
forming on either side of the droplet. These dendrites grow along the droplet-vapor surface until they have 
completely surrounded the liquid Ga subsequently preventing further crystallization. 

This behavior suggests the presence of a Mullins-Sekerka (MS) instability during crystallization, which 
implies that only perturbations of sufficiently large wavelength experience unstable growth. Therefore, 
droplets must be sufficiently large to accommodate perturbations of appropriate wavelength in order for 
the instability to manifest. The effect of temperature. As flux and droplet radius have on the presence of the 
MS instability is given in Section O but the droplet size dependence may be utilized in order to confirm the 
existence of liquid cores by MS instability experimentally. 

It is important to note that simulations suggest that, unlike those formed by surface nucleation, the 
core-shell structure formed due to the MS instability cannot readily be annealed and recrystallized into an 
epitaxial quantum dot. In this case, shell is in registry with respect to the substrate. Therefore, there are 
no defects along which liquid Ga can move so that they may crystallize upon exposure to the As vapor. 
Crystallization of such liquid cores may still occur, however, if the shell is thin enough to dissolve upon an- 
nealing thus exposing the liquid Ga to the As flux, however instability-induced shells are characteristically 
thicker than those resulting from surface nucleation according to simulations. 

Experimental Validation. In order to experimentally validate the presence of the MS instability predicted 
by simulations, nano-islands were formed via droplet epitaxy by depositing Ga droplets of extremely dif- 
ferent sizes (60 and 250 nm diameter) and crystallizing them at high As BEP at 150 °C. Simulations predict 
a size driven transition from a uniform crystallization for small droplet radii to a Ga rich core at higher 
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Figure 10: Snapshots of liquid core formation at times t = 1.4,1.8 and 2.0 seconds after crystallization at 
temperature T = 350°C and deposition rate f^g = 1.0 ML/ sec. 

droplet sizes. We used Si as the growth substrate. This is dictated by the need of differentiation of the nano- 
island constituents (Ga and As) from that of the substrate (Si) [16J, thus allowing the use of standard EDS 
STEM techniques for the characterization of the As and Ga distribution within the islands. The samples 
were prepared as follows. Before the introduction into the MBE system, Si(OOl) substrates were cleaned 
by standard RCA treatment and finally dipped into HE solution to get the H-termination of the surface, 
confirmed by the (1 x 1) pattern observed with RHEED. Subsequently, the substrate temperature was set 
at 780 °C and the hydrogen desorption was carried out until a mixed (2 x 1) (1 x 2) surface reconstruction 
was clearly obtained. After this step, the substrate temperature was decreased to either 200 °C (sample A) 
or 600 °C (sample B) for Ga deposition. Here a Ga molecular beam flux was supplied with a deposition 
rate of 0.075 ML/ s and a background pressure below 5 x 10~^^ Torr for a total of 3.0 ML of Ga. Einally 
the crystallization of Ga was achieved at 150 °C by exposure to an As flux of 5 x 10~^ Torr for 5 minutes. 
During the As irradiation, the RHEED pattern turned from halo, indicative of liquid droplets, to spotty, 
signalling the formation of three dimensional structures. However, a clear spotty pattern was not observed 
for the larger islands. 

The TEM top views of islands from sample A, obtained with diffraction vector g = [220] clearly show 
Moire fringes, due to the interference between the crystal lattice of the GaAs island and that the Si substrate 
(Figure [TTJ. EDS map scans show the same distribution of the Ga and As signals, which mimics the island 
shape, thus demonstrating the correct stoichiometry of the GaAs all over the island volume - i.e. a fully 
crystallized island. A different scenario is shown by sample B. Here the nanostructure generated Moire 
interference fringes are visible only at the island edges with a featureless center. EDS maps scans show 
that the distribution of Ga and As are extremely different inside the island (Figure [TJ]). While the Ga signal 
intensity follows the island profile. As intensity is peaked at the island perimeter with a small, though non- 
zero. As signal coming from the center. Such a distribution is confirmed by EF-TEM (energy filtered TEM) 
measurements, and implies a metallic Ga core similar to the one depicted in Figure [TOl 

The experiment rules out the presence of a liquid Ga core/ GaAs shell structure by means of surface 
nucleation due to the fact that the shell is monocrystalline and appears only for the high droplet radius. 
Nucleation at the liquid-solid interface is predicted to be driven more by the growth temperature and As 
flux. Because both samples were performed at the same growth conditions, if nucleation driven GaAs shells 
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Figure 11: Sample A. a) Planar view TEM micrograph of two small (70 nm diameter) GaAs islands, b) and 
c) EDS maps of Ga and As, respectively, of the islands in (a). 
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Figure 12: Sample B. a) SEM image of alarge (250 nm diameter) GaAs island. b)Planar view TEM micro- 
graph of a GaAs microisland. c) and d)EDS maps of Ga and As, respectively,of the island in (b). e) and f) 
EF-TEM image for Ga and As, respectively, of another GaAs microisland (TEM image not shown). 



were to occur, it must occur in both samples. Because this is not the case, we may conclude that the liquid 
cores observed in the experiments are driven by an MS instability. 



5 Model Analysis 

Both simulations and experiments have shown that droplet epitaxy can result in a wide range of morpholo- 
gies depending on the growth conditions. By varying and T, we have established their effect on the 
resulting nanostructures observed in simulations. From this data, it is clear that compact islands form at 
low temperature and, as the As overpressure is reduced or the temperature is increased the nanostructures 
become rings. In cases of large f ^s, core-shell structures are observed. The simulation results are summa- 
rized in the structural map given in Figure [l3l In this section, we will appeal to physical and mathematical 
arguments to further explain the simulation results. In particular, the solid lines in Figure [13] that delineate 
the morphological structure will be derived in this section. 

We argue that various morphologies observed both in experiments and simulations can be explained in 
the context of three key processes active during nanostructure formation. In the first Ga atoms in the liquid 
drop will be "wicked'' out of the droplet onto the substrate by capillary-type forces when exposed to an As 
overpressure. These forces arise as it is energetically favorable for As atoms on the surface to become fully 
coordinated with Ga atoms. Next, As atoms deposited near or on the droplet diffuse rapidly (T8| through 
the liquid and attach on a growing GaAs front at the liquid-solid interface, crystallizing the droplet epitax- 
ially. Finally, As atoms may also nucleate near the vapor-liquid interface. These three processes: wicking, 
crystallization and nucleation are illustrated in Figure [lH The relative rates at which these processes occur 
depend on growth conditions and will determine the resulting morphology. 
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Figure 13: Morphological Dependence on Growth Conditions. This nanostructural phase map summarizes 
simulation results of droplet epitaxy and crystallization at various As deposition rates and temperatures. 



5.1 Nanoring Formation 

The formation of nanorings is a competition between the wicking and crystallization processes. The crys- 
tallization process results in GaAs fronts that grow until the liquid Ga is consumed. If the fronts coalesce 
before this occurs, the resulting nanostructure is a compact quantum dot. If however, the liquid Ga is 
consumed before the fronts come together, nanorings result. 

To compare the rates of the two processes, we establish expressions for the velocity of the GaAs front 
as well as the rate at which Ga atoms are wicked out of the droplet. Assuming quasi-static deposition, the 
velocity Vn of the GaAs front under an As deposition rate f scales according to that rate: 



where i is the atomic spacing of the lattice and g is factor that depends shape of the liquid region. In the 
wicking process, Ga atoms are driven from the droplet by the deposition of As on the surface away from 
the droplet. The time scale of the wicking process is given by 




Nucleation 



Figure 14: Schematic of kinetic processes that determine nanostructural development. 
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while the diffusion length scale is given by 




(5) 
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where Dq^ is the diffusion coefficient of Ga diffusing on a mostly Ga terminated surface, in the presence of 
As deposition. We have 



-Rq exp 



kBT 



where Rq = 10^^ sec~^ was defined in Section |2] to be a rate pref actor. The energy barrier Eq^ = 1.26 eV 
describing the diffusion of Ga on surface not purely Ga terminated is obtained from the harmonic average 
of the diffusion rates of Ga-on-Ga diffusion and Ga-on-As diffusion. 

Equations dU and ^ yield a velocity v^^ for the wicking of Ga atom from the liquid droplet 



A 



W^oFas exp 



F' 



Ga 



2k bT 



(6) 



A fully crystallized quantum dot forms under conditions where the crystallization process is dominant 
(vn ^ Via). However, if the wicking process is sufficiently fast {vn ^ Vw), then as described above, the fronts 
may fail to coalesce, resulting in a nanoring. The critical configuration separating these two scenarios occurs 
when the crystallization fronts are tangent to each other. To express this critical condition quantitatively, 
consider the volume of unconsumed liquid Ga at time t\ 



V{t) = VQ-nv^t-lpi{gVntf 



(7) 



where pi = ^ is a shape constant describing the geometry of the crystallization front (which we model as a 
quarter-circle) and 

is the initial volume of the liquid droplet of radius r, p2 = ^ being the shape constant that describes the 
geometry of the droplet. As a first-order approximation, the value of p2 is selected as the average of the 
shape constant describing a circular and half -circular droplet. The radius of a Ga droplet obtained from the 
deposition of 4.0 ML of Ga atoms at temperature T and flux Fq^ was found empirically from independent 
simulations to follow the power law: 

where Tq = 11.34 (# atoms) and oc = 0.182 are empirically obtained values. Dq^ is the diffusion coefficient 
for Ga-on-Ga diffusion with corresponding energy barrier Ega = 0-9 eV. The positive root if oi equation 
(O describes the time when all Ga has been consumed. The length of the crystallized front is then given by 
Vntf, and the critical condition may be expressed as 

Vntf = r 

Using equations and (O , this critical configuration can be written in terms of the growth conditions 

as 



{p2-2g^pi)ro 



exp 



Ga 



(9) 



If F^s > F^g then compact structures will result on other hand if F^g < F^^ nanorings will result. Figure 
[Zl illustrates the transition from nanorings to compact quantum dots with increasing F^g. The black line in 
Figure [131 shows a plot of F^^ vs. 1/T. The agreement with the simulation data is quite good - it separates 
compact structures (black dots in the figure) from nanorings (black circles). 

A crucial assumption has been made in the argument presented here, namely that no nucleation takes 
place within the Ga droplet and the crystallizing front moves in a stable fashion (i.e. no Mullins-Sekerka 
instability occurs). In what follows we shall examine in what regimes in parameter space these factors play 
an important role. 
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5.2 Nucleation 



Here we will demonstrate that when the As deposition rate is sufficiently large and temperature is low 
enough, the dominant process will be nucleation at the vapor-liquid interface. As a consequence, the wick- 
ing process makes a negligible contribution to the morphology. The rate of the nucleation process at the 
vapor-liquid interface may be estimated by considering the As concentration c^s(x,y, t) within a domain 
of liquid Ga in contact with a GaAs substrate and As flux. For simplicity we will consider the domain 
to be rectangular with periodic boundary conditions in the horizontal direction. Assuming quasi-static 
deposition at a rate f ML/ sec and temperature T, the concentration satisfies 
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dCAs 
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CAs 



y=0 



0, 

pAsf 

Co 



where Cq is the equilibrium As concentration above a flat liquid-solid interface and the diffusion coefficient 

D As = (^Roexp 



describes the diffusion of As through liquid Ga, = 0.7 eV being the energy barrier for diffusion through- 
out the liquid droplet defined earlier in Section |2l 
The equilibrium concentration Cq is of the form 



Co 



exp 



~ ^A 



where Ejj, are the energy barriers for the detachment and attachment of an As atom in the liquid phase 
from and onto a flat interface, respectively, defined in Section |2l From this diffusion model, the equilibrium 
concentration c^^^^ at the surface of the droplet is given by 
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+ e ^exp 



Ed 



(10) 



where L is the height of the droplet. 

Nucleation is most likely to occur where the As concentration is the largest, near the vapor-liquid inter- 
face. This means that nucleation will occur when Cg^^-j^ is larger than some critical concentration: 



Csurf > Co exp 



(11) 



where Enuc = 0.01 eV is the nucleation barrier of GaAs in liquid Ga, treated as a fitting parameter. Replacing 
L in equation (p^ with the droplet radius r, the critical condition ((TTJ for surface nucleation can be expressed 
in terms of T and F^s as: 
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(12) 



If the As deposition rate exceeds f^^ then GaAs crystallites will form at the vapor-liquid interface, as il- 
lustrated in Figure [HI The blue line in Figure [131 is a plot f^^ as function of 1/T. It accurately predicts the 
presence of nucleation in simulation results (blue triangles). 
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5.3 Mullins-Sekerka Instability 

In the case where the rates of both the wicking and nucleation processes are negligible, the crystallization 
process dominates, resulting in the growth of crystallization fronts at the triple point and in registry with 



As observed in Section 14.2.21 this growth can be unstable due to a Mullins-Sekerka instability leading 
to GaAs shells epitaxial to the substrate surrounding a liquid Ga core. Perturbations to the planar growth 
front of sufficiently long wavelength experience this instability, and a standard linear perturbation analysis 
I34[|35l yields a critical wavelength 



where 7 = 0.1 eV is the liquid/solid interfacial energy, obtained directly from the model. In order to 
accommodate perturbations that experience the MS instability, droplets necessarily must have radius on 
the order of A^, i.e. R > CA^, for some constant C. This critical condition may be expressed in terms of F^g 
and T using the above equation along with the model for droplet radius ©I 



If f ^1 < f < f then liquid cores will form via a Mullins-Sekerka instability. Figure [13] shows a plot 
of f ^1 vs 1 / T (with C = ^ ) . The simulation results when a liquid core was observed are plotted as red 
triangles. We observe that the theoretical curve slightly underestimates the instability within the simula- 
tions. This underestimation is inherent to our model and can be attributed to discrete effects. Such effects 
on nucleation and hence instabilities of the type outlined above are indeed well-studied in the context of 
diffusion limited aggregation 113611371 . 

6 Conclusion 

We have presented experimental, simulation and analytical results detailing the precise relationship be- 
tween growth conditions and the resulting nanostructures formed in droplet epitaxy and crystallization 
experiments. The KMC model used in conjunction with experiments was presented as an explicitly atom- 
istic, multi-species, multi-phase model capable of simulating all the relevant processes involved with GaAs 
homoepitaxy, droplet epitaxy and crystallization. As an initial validation, the KMC model accurately re- 
produces surface termination diagram in the case of GaAs homoepitaxy. We then presented simulation and 
experimental results of droplet formation and crystallization, exhibiting a qualitative agreement on the re- 
sulting morphological dependence on As flux and temperature. Both experiments and simulation suggest 
a continuum of structures ranging between compact quantum dots to broad nanorings as a function of As 
flux and temperature. Simulations also suggest the existence of Ga/GaAs core shell structures, as well as 
elucidating the mechanisms behind their formation. The simulations suggested that a quantum dot with 
a liquid core could result due to a Mullins-Sekerka instability during the growth of the GaAs front. We 
presented experiments which validate this hypothesis. In the case of nucleation-driven shells, we provide 
simulation evidence suggesting that under high temperature annealing, polycrystalline shells may be re- 
crystallized into an epitaxial GaAs island. Lastly, we developed a unifying theory identifying three key 
processes active during crystallization. It predicts the existence of all the phenomena observed above and 
their dependence on growth conditions and compares well with empirical simulation data. 
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